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We attempt to solve the magnetic structure of the gadolinium analogue of 'spin-ice', using a mix- 
ture of experimental and theoretical assumptions. The eventual predictions are essentially consistent 
with both the Mossbauer and neutron measurements but are unrelated to previous proposals. We 
find two possible distinct states, one of which is coplanar and the other is fully three-dimensional. 
We predict that close to the initial transition the preferred state is coplanar but that at the low- 
est temperature the ground-state becomes fully three-dimensional. Unfortunately the energetics 
are consequently complicated. There is a dominant nearest-neighbour Heisenberg interaction but 
then a compromise solution for lifting the final degeneracy resulting from a competition between 
longer-range Heisenberg interactions and direct dipolar interactions on similar energy scales. 

PACS numbers: 03.75.Lm, 39.25.+k, 67.40.-w 



I. INTRODUCTION 

'Spin-ice' is now a well studied experimental system 
with interesting fundamental aspectspj. Originally the 
experimental investigation was undertaken to try to shed 
light on magnetic frustration. The lattice geometry is 
that of a pyrochlore: a three-dimensional network of 
corner-sharing tetrahedra. Although the triangle is ge- 
ometrically frustrated a tetrahedron is much more so, 
culminating in an extensive degeneracy for the antiferro- 
magnetic ground-state manifold for a Heisenberg bonded 
pyrochlore Q. 'Spin-ice' has the required geometry, but 
has a strong orbital and dipolar character that mean that 
the interactions are not Heisenberg-like and the ground- 
state is not even locally antiferromagnetic. This 'failure' 
has led to the current interest [![, but the original issue of 
pyrochlore magnetism remained. Other rare-earth atoms 
are much more Heisenberg-like, with gadolinium being 
archetypal. It is spherically symmetric as a consequence 
of having a half filled shell and, therefore, can be modeled 
as just a large fairly classical spin. The analogue 'spin- 
ice' material Gd2Ti207[3j has now been studied and the 
original issue of a Heisenberg pyrochlore has been inves- 
tigated, subject to any residual dipolar issues. 

The magnetic structure determination for Gd2Ti207 
and Gd2Sn207 is currently a mess. The neutron scat- 
terers predict bizarre states [HHI, with different spins or- 
dering with different moments and the system making 
no attempt to minimise the natural interactions. In con- 
trast, the Mossbauer studies |6[ offer a simple picture of 
all sites equally ordered and the additional restriction 
that the moments lie in a particular plane perpendicular 
to a natural local crystallographic direction. In this ar- 
ticle we try to rationalise both types of experiment and 
achieve a consistent magnetic structure. 

The neutron scattering provides a magnetic state in- 
dexed by 5j 5^ • Since the system has an intrinsic four 

atoms per unit cell, this leads to thirty-two atoms per 
magnetic cell, rather more than is comfortable. We con- 
sequently start from the Mossbauer experiments which 
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FIG. 1: Restriction of spins to a plane perpendicular to the 
local crystallographic directions pointing to the centre of the 
tetrahedron. 



suggest that at low temperatures all moments are ordered 
with the restriction that they are ordered perpendicular 
to the natural local crystallographic direction. (Note that 
for spin-ice the spins are parallel to these local crystallo- 
graphic directions, which corresponds to a reversal of the 
sign of the anisotropy interaction.) For a single tetra- 
hedron this restriction is depicted in Fig[TJ We stress 
that this restriction is not expected to be a dominant en- 
ergetic restriction because the spin-orbit/dipolar aspects 
are likely to be very weak in comparison to the Heisen- 
berg interactions. We are employing these constraints as 
experimental restrictions which only become fulfilled to 
extract some extra degeneracy from the choice of ground- 
state, once all stronger Heisenberg interactions have been 
optimised. Physically our arguments are consequently 
dangerous since slight avoidance of our constraints can 
provide energetic savings but be lost in the noise of the 
experiments! 
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We further enforce an energetic constraint that the to- 
tal spin in each tetrahedron vanishes, which is suggested 
by the experiment but is not guaranteed by it. This then 
provides a tractable problem, which annotates the resid- 
ual degeneracy of the Heisenberg model subject to the 
additional experimentally observed single-ion anisotropy. 
Note that this antiferromagnetic constraint is not appli- 
cable to spin-ice[l|, where the orbital and dipolar com- 
plications provide a very non-Heisenberg interaction. 



II. MAGNETIC SOLUTIONS 




The solution of our constrained problem is much sim- 
plified using the normalised basis: 
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[x cos Xq + y cos yo + z cos zq] (1) 



[xcosici — y cosj/i — zcoszi] (2) 




FIG. 2: The six possible angles for a particular global spin 
configuration. This example is So and each spin configuration 
is described simultaneously by all three angles, xq, yo and zq. 



-X COS X2 + y COS J/2 — Z COS Z2\ 
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and then three associated with pairs 



all modulo 2-7r. This guarantees that the spins lie in the 
appropriate planes and also that they are normalised. 
The spins may also be decomposed as: 
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with analogues for the other three spins, which relates 
the chosen angles to the underlying Cartesian basis. 
The solutions to the constraint (see Appendix) 



S + Si+S 2 + S 3 =0 



(9) 



may be characterised by the different ways of pairing the 
atoms. We have one style of solutions with all the angles 
equal 

Xq = Xi = X 2 =X 3 2/0=2/1= 2/2 = 2/3 

Zq = Z\ = Z 2 = Z Z (10) 



Xq = X\ = —X-i 
ZQ = -Zl = 



-x 3 2/o 

-Z2 = Z 3 



-2/1 = 2/2 



-2/3 



fill 



for the three choices of pairs of pairs. 

The restrictions to the anisotropy plane and to antifcr- 
romagnetism severely constrain the permitted magnetic 
states. Once we have fixed one spin then there are (a 
maximum of) six permitted orientations for each spin 
in the local basis (a maximum of twenty- four globally). 
These states are generated by x a h-> —x a , y a h-> —y a and 
z a i-> — z a , which are the only possibilities permitted by 
the Heisenberg model. The possible states are annotated 
in Fig[2] One can then generate possible global states 
from the initial spin, using one out of the four possi- 
ble configurations in each subsequent tetrahedron, as de- 
picted in Fig|31 together with the analogues with n <— > n 
and independently the inverted tetrahedra. 

The next step is to enumerate the possible global con- 
figurations. This proves quite difficult, however, so first 
we set the scene. For spin-ice the spin orientations are 
quite different, being controlled by 'two in and two out' 
for each tetrahedron. The spins are oriented towards 
the centres of the tetrahedra which induces less technical 
problems. All possible ground-states may be generated 
by propagation parallel to a Cartesian direction. One 
can place an arbitrary collection of spins in a particular 
plane of the geometry which places two spins in each par- 
ticipating tetrahedron. Any pairs that are either both in 
or both out propagate uniquely and each pair that are 
one in and one out provide two possible continuations. 
For each possible state we then propagate to the next 
plane generating an exponentially large number of pos- 
sible ground-states. Note that this is not a solution to 
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FIG. 3: Possible spin configurations in an optimised tetrahe- 
dron, using the notation of FigO 



the spin-ice-like anisotropy projected Heisenberg model, 
which finds all spins in or all spins out and hence only 
two distinct ground-states because when one spin is cho- 
sen then all others are constrained. 

For our model we follow a two stage process. Firstly we 
ignore the 'bar' and look only at the integer n € {1, 2, 3}. 
For the special case of Zo=0 this distinction becomes 
moot anyway. Indeed, if we think a bit more physi- 
cally, then we would not expect the spins to employ arbi- 
trary directions and the special cases of z e {0, 4^} 
and Zq £ { — §,§,§} offer distinct extra symmetries with 
both involving only twelve global spin orientations. The 
first provides a restriction to only three local orientations 
and the second finds commonality between distinct sub- 
lattice configurations leading to the orientations: 



±y ± z ±z ± x ±x ± y 

l To ' To ' To >' 
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(12) 



This second option is probably the best candidate for the 
experimental systems. By employing propagation along 
one of the Cartesian directions, we can use a similar ar- 
gument to the previous to investigate the degeneracy. 
We focus on a plane perpendicular to the z-direction for 
which Fig|3]provides the restrictions that 1 and 2 may not 
be neighbours in the plane but pairs of 1 and 2 may prop- 
agate to the next plane using either 1 or 2 for the next 
pair. This provides a large but seemingly not exponential 
degeneracy, since as well as the new configurations from 
the choices of 1 and 2, there is also an associated loss from 
discounting states where 1 and 2 are neighbours in the 
next plane. The inclusion of the bar further complicates 
matters. 

We now consider the scattering experiments and in- 
clude the idea that the scattering can be indexed using 

( 2 ' 5 ' 5 ) • ^ tne simplest level this restriction amounts 













2 






3 








5 













i 






4 














1 






7 













14 



FIG. 4: Distinct atoms in the magnetic unit cell, n are in the 
reverse direction to n. 
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FIG. 5: Possible spin configurations for the case «o=-S- 
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to choosing the sixteen spins of FigHJ subject to the pre- 
vious restrictions of Fig(3J combined with the constraint 
that the 'barred' spins are reversed in orientation. In 
our modeling, this restriction to having reversed spins is 
also a strong constraint. The easiest way to enforce it is 
to restrict attention to ^o=f and then in Fig J5] the two 
spins n and n are actually in opposite physical directions. 
The previous problem of annotating permissible configu- 
rations in terms of the states n now provides an elemen- 
tary solution: we find four styles. Firstly, pure solutions 
where each site has the same n. Secondly, alternating 
solutions where in a particular direction planes alternate 
between two distinct values of n. Thirdly, tetrahedral so- 
lutions with alternating tetrahedra alternating between 
two distinct values of n. Fourthly, period-four solutions 
where one plane in four has a distinct value of n. Figj5] 
shows the only two examples which are are consistent 

with 

A second set of solutions may be obtained from z =-^ 
where now the pairs {1,2} {2,3} and {3,1} correspond 
to reversed spins. The permitted solutions are depicted 
in FigjBJ and there are only four physically distinct styles. 
Clearly by employing symmetry there are a mass of ro- 
tation and translationally associated analogues to the 
depicted states. In total, however, there appear to be 
only six styles of solutions. Careful investigation of these 
states demonstrate that in fact there are only two classes 
unrelated by symmetry, which we can choose to be the 
zo=f states of FigJS] The first is coplanar and the second 
fully three-dimensional. 
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FIG. 6: Possible spin configurations for the case zq—-j. 



FIG. 7: Kagome plane decomposed into the two inversion re- 
lated triangles. The (a) denote atoms that complete tetrahe- 
dra above and the (b) denote atoms that complete tetrahedra 
below. 



III. EXPERIMENTAL AND THEORETICAL 
ISSUES 

There is one critical final piece in the puzzle: the form- 
factors of the magnetic Bragg spots. The internal struc- 
ture of our proposed states can be probed by measuring 
the relative intensities of the different magnetic Bragg 
spots. In our notation the strengths of the nearest Bragg 
spots to the origin are directly related to the magnetic 
moments on each of the four underlying sublattices: 

Bo = S + S 1 + S 2 + S 3 Bi = S 4 + S 5 + S 6 + S 7 (13) 

B2 = Sg + Sg + S10 + Six B3 = S12 + S13 + S14 + S15. 

In all of our states these quantities vanish! 

The actual Bragg spot, and its symmetrically 

related neighbours, correspond to decomposing the lat- 
tice into alternating planes of Kagome connectivity and 
sparse triangular planes, keeping the phase within such a 
plane uniform and alternating the phase between neigh- 
bouring equivalent geometries. Interestingly, our antifcr- 
romagnetic ansatz directly controls this character. The 
Kagome planes may be decomposed into two inversion 
related sets of triangles, as depicted in Fig{7] One set 
(a) corresponds to a set of triangular faces of tetrahedra 
which are completed above the plane, and the other (b) 
form a set of corresponding tetrahedra which are com- 
pleted below the plane. These completing atoms form 
the neighbouring sparse triangular planes. The imposed 
constraint that each tetrahedron has zero total spin then 
enforces that the total spin in a Kagome plane is an- 
tiparallel to both neighbouring triangular planes. The 
phase relationship that 'neighbouring' triangular planes 



should be antiparallel then requires that the total spin in 
each Kagome plane must vanish. Consequently the asso- 
ciated Bragg spots must also vanish. Physically the argu- 
ment only depends on the nearest-neighbour Heisenberg 
interactions: solving each tertrahedron gives the require- 
ment that that tetrahedron should have vanishing total 
spin. The Kagome argument of FigjT] requires that the 
total spin of each neighbouring plane along the (1,1,1) 
direction (or analogues) must be antiparallel. Since the 
relevant Bragg spots add every second plane with an al- 
ternating sign, the appearance of such a Bragg spot is 
inconsistent with the nearest-neighbour Heisenberg inter- 
action. We have also imposed the additional constraint 
that our solutions have the observed magnetic periodicity 
which requires that the total-spin of each plane vanishes. 

This argument provides a second important physical 
consequence. In terms of the original four atoms per unit 
cell of the underlying pyrochlore structure, the triangular 
lattices between Kagome planes provide each of the four 
sublattices as the orientation of the Kagome planes is var- 
ied. This fact tells us that the total spin of each sublat- 
tice independently vanishes. Although we have developed 
our argument employing spin anisotropy and experimen- 
tal periodicity, our chosen states automatically minimise 
spin interactions that we did not impose. To see this we 
next analyse third-nearest neighbour Heisenberg interac- 
tions, which in this case are expected to be stronger than 
the second-nearest neighbour interactions due to the na- 
ture of relavent interaction pathways. We find that the 
third-nearest neighbour interactions are restricted to con- 
nect atoms on the same original sublattice. There are 
twelve such third-neighbours to each spin and they come 
in two classes. Firstly, there are six neighbours with a 
magnetic gadolinium atom lying exactly in the middle 
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FIG. 8: Cage of gadolinium atoms surrounding a tetrahedron 
of titanium atoms. 

between them. Secondly, the other six neighbours have 
a non-magnetic titanium atom lying between them. The 
gadolinium chains form the original face-centre-cubic lat- 
tice but with bonds omitted between neighbouring atoms 
lying in planes perpendicular to the appropriate (1,1,1) 
direction. In our notation, this interaction is optimised 
by maximising the total spin of the four atoms on the 
same sublattice rather than minimising it. The chains 
of alternating gadolinium and titanium atoms yield the 

I 



r k = 

with, 

I/* l/i If* 

x = cos — y = cos z = cos (15) 

2 2 2 v ' 

ky + k z k z + K x k x + ky 

X = cos — - Y = cos Z = cos - 

2 2 2 

(16) 

Ci=COsfc x C2 = COS ky C3 = COSfc 2 (17) 

I 



sparse triangular lattices in-between the Kagome planes. 
Although our vanishing total spin on these planes is a low 
energy state for these third neighbour interactions, it is 
only the ground-state of the triangular lattice when there 
are additional longer-neighbour interactions also present. 

We now embark upon a low level analysis of the plau- 
sible energetic interactions responsible for our proposed 
states. There is a sizable literature which we do not do 
justice to. The most relevant mean-held analysis of the 
Heisenberg interactions is not compatible with the actual 
interactions in this systemQ, using physical distance in- 
stead of active pathways to choose the relevant interac- 
tions. The dipolar investigations are more relevant [1, ||| 
and indicate that the lowest energy solution should be a 
q=0 state 3 although (q,q,q) states should also be low 
energy states^ and therefore competitive. We employ 
much simpler arguments to try to get to grips with the 
most likely physical explanation for the observed experi- 
mental states. 

The actual additional Heisenberg interactions proba- 
bly stem from pathways across titanium atoms, since the 
hopping matrix elements are so much larger. The cage 
of atoms depicted in FigJS] are all expected to develop 
Heisenberg interactions of various strengths between the 
members, although the bonds between members of the 
same sublattice should be largest. We employ two new 
matrix elements for the additional Heisenberg interac- 
tions across the hexagons depicted in Figj8] A for dia- 
metric and k\ for the remaining two interactions. The 
structure factor is then controlled by the matrix: 



(14) 



I 

as the associated parameterisation. In the absence of the 
perturbations this structure factor has two degenerate 
bands at T — — 1 which control the spin degeneracy. We 
treat the additional perturbations as small and solve for 
the lifting of this degeneracy in the two degenerate bands. 
Employing the parameter e to control the eigenvalues: 

r = -l-A-2Ae, (18) 

we hnd that the perturbed structure factor satisfies: 



"A (4xyz - 1) X + 2Xkxci 

X + 2Xnxc x A {AxYZ - 1) 

Y + 2\nyC2 z + 2\kZcz 

_ Z + 2Xkzc3 y + 2XkYc2 



Y + 2Xnyc2 Z + 2Xkzc3 

z + 2XkZcs y + 2XkYc2 

X (4:XyZ — 1) x + 2XnX Cl 

x + 2XkXc x X(AXYz-1)_ 



(l-a 2 )e 2 -2{a2[l+AK]+?>al-a l [Za 1 +a3][l+K])e+Za 2 ^^^ = 

(19) 
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FIG. 9: Angle of the spiral measured in units of n. 



in terms of 



ax 



ci + c 2 + c 3 



C 2 C 3 + C3C1 + CiC 2 



02 = ^tt z — a 3 = C1C2C3. 

6 3 

(20) 

This can then be minimised to predict the likely posi- 
tions of the Bragg spots stemming from these physical 
interactions. Although the observed experimental Bragg 
scattering is at low energy, there is always a lower energy 
solution that finds a spiral parallel to one of the Carte- 
sian axes. The ground-state as a function of k involves a 
solution of the form: 



ci = c 2 



cos 8 C3 = 1 . 



(21) 



This spiral starts out with 9 = ^ when k — and ends up 
at the experimental 8 = ^ with the unphysical value of 
K = 1, as depicted is in FigJHl The energy of this ground- 
state is compared to that of the experimental solutions 
and the zone-centre solution in FigfTOl It is clear that the 
second-neighbour Heisenberg interactions are not suffi- 
cient to describe the physics, and that the zone-centre 
state has a very poor energy. It therefore appears neces- 
sary to resort to the dipolar interactions to energetically 
explain the observed phases. 

Dipolar interactions take the generic form 



H = 



Sq-Sff — 3S a .r a pr a p.S/3 

M~ 3 



(22) 



where r a p = r a — rp is the vector connecting the 
two interacting spins. This interaction is mathemat- 
ically taxing because it is both long-range and it re- 
lates the spin orientation to the lattice directions which 
breaks spin isotropy. We will extract the initial isotropic 
Heisenberg interactions which renormalise the existing 
exchange-based interactions and focus only on the second 



FIG. 10: Energies of various states, ground-state (bold), 
experimentally observed state (dotted) and zone-centre 
(dashed) . 
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FIG. 11: Total-spin zero eigenstates of the short-range dipolar 
model, in order of energy. 



term. This second term we further restrict to nearest- 
neighbours: 



H 



-3ES 

a>/3 



(23) 



It is elementary to solve this rescaled dipolar interaction 
on a single tetrahedron and we find five styles of solution 
at energies: {-|(1 + ^17), -3, 3, §(^17 - 1), 12}. The 
ground-state has 86.4% of its moments parallel and is 
consequently a high energy state if the nearest-neighbour 
Heisenberg interaction dominates. The next state is 
triply degenerate and as it has zero total-spin it is com- 
patible with the Heisenberg interaction. The third state 
is doubly degenerate and also has zero total-spin. The 
fourth state is not compatible as it has 13.6% of its mo- 
ments parallel. Finally the fifth state also has zero total- 
spin. The three compatible styles of states are depicted 
in Fig QT] and the crucial observation is that the two low 
energy states involve spins which are perpendicular to 
the natural crystallographic directions (i.e. the spins are 
in the planes indicated by the Mossbauer experiments). 

The planar spin restrictions deduced from the 
Mossbauer experiments can be derived from the nearest- 
neighbour dipolar interactions subjected to dominant 
nearest-neighbour Heisenberg interactions. Enforcing 



7 



equal length spins and including the Heisenberg con- 
straints, Eq.@, we may rewrite the dipolar interactions, 
Eq.((23j), as 

H = l [(S$ + S» + Stf + {St - Si - S{) 2 
+(-Sf + SI - S z 2 f + (-Sf - SI + Stf] 

+(sy~s! + sy-syf 

+ (5 2 - St - S z 2 + SI) 2 ] . (24) 

This is clearly minimised by zeroing each of the seven 
quadratics, which is consistent and indeed uniquely spec- 
ifies the three states shown later in Fig. In particular 
we note that the first four quadratics may be represented 
by 

[S .(l,l,l))] 2 [Si.(l,-l,-l)] a 
[S 2 .(-l,l,-l)f [S 3 .(-l,-l,l)] a • (25) 

Orienting the four spins perpendicular to their local crys- 
tallographic directions therefore partially minimises the 
local dipolar energy. 

One might now conclude that the problem is solved. 
We have a dominant nearest-neighbour Heisenberg in- 
teraction with a weaker dipolar interaction which sta- 
bilises the observed experimental states, but this is not 
the case. There is a clear solution to the problem of dom- 
inant nearest-neighbour Heisenberg and weak nearest- 
neighbour dipolar forces and it is not the experimentally 
observed state. Instead all that is required is to repeat 
the spiral solution of Fig llll in all directions to provide a 
compatible solution that has k = 0@. The only way to 
now rationalise the experiments is to assume that there is 
a competition between the dipolar forces and the longer- 
range Heisenberg interactions. The energy of the Heisen- 
berg interactions at k = is e=-2+n which is much worse 
than e=0 for the experimental solutions. We can then 
argue that dipolar interactions are not strong enough to 
overcome this energy but are strong enough to select the 
observed ground-state over the weak preference by the 
second and third-neighbour Heisenberg interactions. 

The next task is to assess the dipolar energies in the 
language that we have used to describe our solutions. 
The dipolar states which do not respect zero total-spin 
can be split into pieces which are parallel and perpendic- 
ular to the natural crystallographic directions. We find 
the 'two in and two out' state of spin-ice parallel to the 
crystallographic directions and a saturated ferromagnetic 
state parallel to a Cartesian direction. These states are 
consequently totally irrelevant. The high energy state 
is also restricted to the crystallographic directions and 
hence we are restricted to the first two states depicted in 
FigCH Careful calculation of the dipolar energies of the 



six bonds in a tetrahedron then yields 

D = cos Xq cos Xl + cos X2 COS £3 

+ cos j/o cos j/2 + cos ?/i cos 2/3 (26) 

+ COS Zq COS Z3 + COS Z\ COS Z2 

and when we further incorporate our solutions to nearest- 
neighbour Heisenberg interactions we find 

D = 3 D\ = 2,cos{2x a ) 

D 2 = 3 cos(2j/ Q ) L> 3 = 3 cos(2z Q ) (27) 

for each of the configurations depicted in the columns of 
FigEl (where a € {0,1,2,3} is any of the spins). Min- 
imising these energies clearly leads to two of the general 
dipolar eigenstates (shown in Fig. [TTj) . However, it also 
highlights the solutions zo S { § , ±§ } that we previously 

required, in order to be consistent with i, ij. We 

can now, therefore, choose the value of zq to optimise the 
dipolar interactions rather than to fit the experiments. 

The final task in this section is to analyse the solutions 
that we previously found for their dipolar content. The 
first half of the states depicted in both FigJS] and Figj6] 
involve half of the tetrahedra being in the dipolar ground- 
state and the other half being in the high energy state. 
The second set of states also involve half of the tetrahe- 
dra being in the dipolar ground-state but the other half 
are not in eigenstates of the dipolar interaction at all, but 
instead are superimposed three-quarters the high energy 
state and a quarter the ground-state. Consequently the 
more exotic states are actually lower in dipolar energy 
and therefore expected to be favoured at low tempera- 
tures. 

One crucial observation is that for our coplanar states 
the dipolar interactions are energetically consistent with 
the nearest-neighbour Heisenberg interactions while for 
the non-coplanar states they are not. Although the spin 
state in each tetrahedron is not always the ground-state, 
for the coplanar state it is always a local eigenstate. For 
the non-coplanar states, however, half of the tetrahedra 
have states which are not local eigenstates. Consequently 
the system would be expected to relax slightly to improve 
the dipolar energy a little at the expense of the nearest- 
neighbour Heisenberg energy. This provides a physical 
mechanism for the magnetic Bragg spots closest to the 
origin to appear! 

IV. STATE CHARACTERISATION 

Our final task is to try to characterise our different 
solutions so that they might be separated both theoret- 
ically and experimentally. The joint issues of quantum 
fluctuations, thermal disorder and static disorder will all 
lift degeneracy and prefer particular states. Quantum 
fluctuations prefer collinear states which allows more 
coherent fluctuations. Static disorder, however, prefers 
non-collinear states [lfj, which allow static distortions 
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in perpendicular directions to the disorder field. Ther- 
mal fluctuations prefer the highest density of states at 
low ener gies , to optimise entropy, so called 'order from 



disorder' |ll|. This tends to offer similar states to those 
stabilised by quantum fluctuations. Our states come in 
two classes, the first half of each set of states in FigEJand 
Fig are coplanar and the second half involve fully three- 
dimensional spin orientations (although admittedly only 
eight of the twelve conceivable). We might anticipate the 
coplanar solutions to be stable at low temperature, on 
quantum fluctuation grounds, but alternatively we might 
expect the non-coplanar states on dipolar grounds. It is 
natural to associate the observed transition in Gd2Ti2 07 
as a transition between the two possible states, but we 
have no overwhelming experimental insights that provide 
clues as to which order of states one might like to pick. 

The experiments do provide an important anomaly. 
Although in the first investigation [4| the Bragg spots 
closest to the origin were seen to vanish, the more care- 
ful investigation !5!j found small intensities which were 
attributed to the low temperature phase. As we have 
pointed out, such spots are inconsistent with the nearest- 
neighbour Heisenberg model. There are various options. 
Firstly, the observed scattering might be diffuse and as- 
sociated with disorder. It is known that diffuse scattering 
is much larger than usual in non-collinear magnets [l2|. 
Secondly, the magnetic state is expected to break the cu- 
bic symmetry and hence there is no longer a requirement 
that the Heisenberg interaction is perfect. Thirdly, the 
scattering might be nuclear in origin and hence be as- 
sociated with a tiny structural distortion (although this 
appears physically quite unlikely). The dipolar forces 
argument suggests that coplanar states yield no Bragg 
intensity but that non-coplanar states yield a small one; 
this is the only clue we have to offer. 

The actual spin arrangements are quite complicated 
and hence are difficult to depict. We have taken half 
of the magnetic cell and pictured the coplanar and non- 
coplanar states respectively in FigH^] and FigO The 
complete spin arrangements are constructed from these 
by A-B ordering a full cubic super-lattice using the de- 
picted cell and its spin inverse as the decoration for the 
super-lattice. 



V. PARADOX 

We now arrive at a theoretical paradox: our calcula- 
tions are inconsistent with the current theoretical litera- 
ture! We have found a complete set of solutions to our 
assumptions and the local dipolar energy is much worse 
than the q = solution that is expected from the 'order 
from disorder' calculations [8|. The entire ethos behind 
that calculation was the lifting of a degeneracy that re- 
mained when the dipolar interaction was included, essen- 
tially exactly in the classical limit at zero temperature. 
In our calculations we ought to have found a solution 
that was degenerate with the q = solution. 




FIG. 12: Coplanar spin arrangement. 




FIG. 13: Non-coplanar spin arrangement. 



There are various possible resolutions. Firstly, we have 
restricted our calculations to states which are orthogonal 
to the natural crystal directions and this eliminates two 
additional zero total-spin states from our optimisation. 
One is the state with spins pointing to the centre of the 
cube and the other is the state depicted in FigfT^J Both of 
these states are of much higher energy than the ground- 
state. Use of these additional states, therefore, cannot lift 
us to the energy of the q = solution. Indeed, we can 
exactly solve the problem of dominant nearest-neighbour 
Heisenberg and local dipolar interactions: we are re- 
stricted to the states depicted in Fig[T5] These states, 
which are in the perpendicular subspace, constitute only 
three of the possible configurations in FigJ3J For the case 
zq = they are the three configurations with {n, n, n, n}, 
and we can only construct q=0 solutions from these con- 
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FIG. 14: Possible spin configuration ignored so far. 
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FIG. 15: Local dipolar ground-states when total-spin of tetra- 
hedron vanishes. 



figurations Secondly, we have only considered the local 
dipolar interaction and the longer-range contributions 
could overturn our argument. Although the dipolar in- 
teraction is long-range, the divergence is irrelevant and 
we do not anticipate issues here. Thirdly there might be 
a loop-hole in the previous theoretical arguments. The 
previous calculations for the Heisenberg and dipolar in- 
teractions combined Q calculated the effective structure 
factor for the classical magnetic problem, but they did 
not proceed on to provide an actual magnetic solution. 
It was presumed that an appropriate magnetic solution 
would exist. Unfortunately, this is not guaranteed. For 
classical magnetism there are additional constraints: the 
magnitude of the spin of each atom on each sublattice 
must independently be normalised (l3j . For Heisenberg- 
like exchange models there is always an elementary spiral 
(in all bar the Ising model) that will provide a solution. 
In this dipolar problem, however, one needs to make a 
multiple-q state involving all four appropriate q's, and 
there are simply not enough spin dimensions to satisfy all 
the constraints. It would appear that solving the dipolar 
problem is harder than first imagined. 



VI. MORE EXPERIMENTAL AND 
THEORETICAL ISSUES 

In this section we range across the main experiments 
that have been performed on Gd2Ti 2 C>7 and comment on 
their relationship to our and previous proposals in the 
literature. We start out with local probes: Mossbauer 
and muon spin rotation. 

The Mossbauer experiments 6] offer the assertion that 
the spins are all the same length and that they are ori- 



ented perpendicular to the local crystallographic direc- 
tions, assumptions that we employed in our modeling. 
More careful reading exhibits the possibility of two mag- 
netic moments not both in the preferred planes, but only 
in the initial intermediate temperature phase. An ad- 
ditional degree of freedom always provides a better fit; 
the improvement is not significant. The second magnetic 
moment is also much too large to agree with the states 
proposed by the neutron scatterers. The low tempera- 
ture phase does not appear to accept a phase with two 
different moments. 

The muon spin- rotation experiment [3] offers two 
clues. Firstly the muon sees two magnetic sites with 
slightly different fields and a non-magnetic site (or a par- 
allel spin). Secondly there is a sizable relaxation even 
at very low temperatures. Due to the lack of knowledge 
about where the muon sits (probably in a low symme- 
try position close to one of the oxygens) the magnetic 
information is not easily usable. The relaxation tells us 
that there are active low energy excitations disturbing 
the spins even at the lowest temperatures. Also present 
in this investigation is some specific heat data, which is 
most instructive. The majority of the entropy is dis- 
sipated above the magnetic transitions and the entropy 
goes to zero at zero temperature with a power law. There 
is no residual degeneracy at zero temperature and the 
low temperature behaviour is well represented by some 
low energy excitations (surprisingly low energy perhaps) . 
It is therefore, not possible that a macroscopic fraction 
of the spins are locally disordered in the vicinity of the 
magnetic phase transition. Any disorder must be highly 
correlated and associated with minimal entropy, viz a few 
well defined excitations. 

There are some very interesting neutron spin-echo 
measurements [l5l Il6j. This pretty experiment gives di- 
rect access to the time dependence of the magnetic cor- 
relations. The long-time limit clearly shows the ex- 
pected long-range order. The temporal decay of the 
measurement exhibits the time-scale on which the fluc- 
tuations disorder the initial 'snapshot' of the spin corre- 
lations. Thermal fluctuations are expected to decay at 
low temperature and so the measurement flattens at low 
temperature[l6j . The most intriguing aspect is the fact 
that on observable time-scales the signal does not satu- 
rate, but appears to converge to three-quarters! This is 
taken to mean that a quarter of the spins are fluctuating, 
but from the specific heat data, they would have to be 
doing this coherently. It is not clear from the papers how 
this normalisation is accomplished, because naively one 
takes a very low temperature and zero field to normalise 
the data, assuming that the system has no residual dy- 
namic fluctuations there. If this result is correct, then it 
is certainly interesting. 

Finally we come to the elastic neutron scattering 0,@, 
with the proposal of states with spins of different average 
length. Although it is stated that the neutron scatter- 
ing promotes only these states, the current paper clearly 
shows that states exist which are consistent with the 
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elastic neutron scattering and the Mossbauer restrictions 
that the spins are essentially of the same length. We now 
consider the theoretical ideas and what they say about 
the proposed states. The states with spins of different 
average lengths require fluctuations to explain the re- 
duction in lengths, either thermal or quantum in nature. 
At low temperature the entropy measurements indicate 
that thermal fluctuations are irrelevant and so we are 
left solely with quantum fluctuations. This is a possible 
explanation, because quantum fluctuations gain energy 
by making spins locally antiparallel in directions perpen- 
dicular to the classical order. They do not require any 
entropy, as they amount to a particular phase for the fluc- 
tuation and not a random one. Low dimensional systems 
are particularly susceptible to quantum fluctuations, as 
are frustrated systems, although larger spins are less sus- 
ceptible as the fractional energy gain is proportional to 
the inverse of the length of the spin. The situation that 
we have is abnormal because we anticipate that there 
are two energy scales: the highest wants all tetrahedra 
to have zero total-spin and the second lifts the residual 
degeneracy and promotes the observed long-range order. 
A study of a single tetrahedron sets the fluctuation scale. 
The classical energy of a single terahedron is —2JS 2 and 
the quantum energy is —2JS 2 — 2JS. This means that 
the quantum fluctuations have access to | of the clas- 
sical energy for our system. Obviously, a fluctuation in 
one tetrahedron corrupts two others and so only a small 
fraction of this energy is available to all tetrahedra simul- 
taneously. In our proposed states we would expect quan- 
tum fluctuations to exist and to reduce the observable 
lengths of the ordered moments. There is no expectation 
that some spins would fluctuate vastly more than others. 
This is a plausible interpretation for the lack of saturation 
in the neutron spin-echo experiments. For the variable- 
length spin states, we need to assume that the quantum 
fluctuations are strong for some spins and weak for oth- 
ers. Once again, this is not impossible, because the idea 
of dimerisation can be extended to that of independent 
tetrahedra with quantum states, at a stretch. The prob- 
lem is that one needs to gain more than one loses. The 
states proposed by the elastic neutron scatterers involve 
some spins with essentially saturated classical moments 
and some with tiny moments. If the moment is saturated 
then it cannot fluctuate. The spins with tiny moments 
are either well separated on sparse triangular planes, or 
combined on particular tetrahedra. The first possibility 
is energetically awful, because no coherent fluctuations 
can develop and incoherent fluctuations gain no energy 
and contradict the specific heat data. The second possi- 
bility is only problematic because the energetics is poor. 
Each fluctuating tetrahedron has access to — 7 J of ex- 
tra energy, but the four connected tetrahedra have three 
correlated saturated spins in the 120° phase and con- 
sequently naturally lose — 24.5J, a very bad deal! The 
variable-length spin phases, therefore, are not consistent 
with the theoretical assumption that the dominant en- 
ergy is the nearest-neighbour Heisenbcrg interaction. 



Experiments on Gd2Sn2 07 indicate that the ground- 
state is actually the dipolar preferred q = phase [1 7] . 
At first sight this is surprising, since the degeneracy is 
on the gadolinium atom and the tin or titanium atom 
appears passive. However, we need to employ higher- 
order exchange paths across this 'passive' atom and if we 
assume that the titanium ci-shell is closer to the chem- 
ical potential than the tin p-shell, then we can clearly 
expect a difference in the longer-range Heisenberg inter- 
actions. This would indicate that for the tin compound 
only nearest-neighbour Heisenberg and dipolar interac- 
tions are required to explain it. 



VII. CONCLUSION 

We have employed a hybrid method to determine plau- 
sible magnetic ground-states for the pyrochlore magnet 
Gd2Tl2 07. We used the theoretical ansatz that the 
nearest-neighbour Heisenberg model is minimised, the 
experimental observations that the magnetic Bragg scat- 
tering is indexed by Qjg'l)' anc ^ that all spins have 
equivalent moments which are oriented perpendicular to 
the local natural crystallographic directions. These as- 
sumptions provide only two distinct classes of solutions 
and we expect the system to exhibit one or both of these 
phases to explain the observed phase diagram. Our best 
guess is that close to the initial transition the thermal 
fluctuations prefer a coplanar state but that at low tem- 
perature the dipolar interactions energetically stabilise a 
non-coplanar state. The observed appearance of a very 
weak intensity to the closest magnetic Bragg spots to the 
origin only in the low temperature phase is consistent 
with this prediction. We await more detailed experimen- 
tal form factor analysis to confirm or deny this proposal. 

The field dependent phase diagram [l^ seems particu- 
larly instructive, but currently is too difficult for us to 
predict. 
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Appendix A 

In this appendix we solve the trigonometric constraints 
cosxo + cosxi — cosa;2 — cosa:3 = (Al) 

cos yo — cos 2/1 + cos 2/2 — cos 2/3 = (A2) 

cos zq — cos z\ — cos z 2 + cos z 3 = 0. (A3) 
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Firstly we break the symmetry and focus on z a , rewriting 
the first two equations provides 



sin 
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(A5) 



and then mixing them offers 



z - z 3 \ . / z + z 3 
sin I : ) sin I — - — 

.., . , zi- z 2 \ I Zi + z 2 
+ v3 sin I — - — I cos 



/o • I Z - Z 3 \ ( Z + Z 3 
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(A6) 



zi - z 2 , 
sin — - — sin 
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-zi + z 2 
2 J "~ V 2 

Including the third original equation as 



, z o — z 3 \ j z + z 3 
cos cos 



0. (A7) 



zi - z 2 
cos l — - — I COS 
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zi + z 2 



(A8) 



2 J \ 2 
we can now eliminate zi and z 2 to provide 

s 2 _ (4 S 2 + - 3) ( [9 - 8s 2 + ] [1 - s 2 _] + 24) = (A9) 
in terms of 

'z ± 23" 



s± = sm 



(A10) 



The final complicated solution is unphysical so we gen- 
erate two independent solutions. Firstly (1) 



zo ~ z 3 







zi - z 2 



(All) 



2 J V 2 

and consequently (modulo 2n) 

zq = z 3 zi = z 2 cos zq = cos zi (A12) 
and the two solutions 

Z = Zi = Z 2 = Z 3 Z = -Zl = ~Z 2 = Z 3 . (A13) 

Secondly (2) 

1 



. 2 z + z 3 
sm 



cos(z + z 3 ) 



(A14) 



and from Eq |A6l and Eq |A7l 



, zi - z 2 \ . ( z + z 3 Zi+ z 2 
sm sm 



0. (A15) 



The solution zi—z 2 provides a subset of the solutions of 
(1), and the new possibility is 



ZQ + Z 3 = Z 1 + Z 2 

leading to two new solutions: (2.i) 

47T 



(A16) 



zo + z 3 — — = Zl + z 2 



> xq + x 3 = = Xi + x 2 
(A17) 

which reduces the original equations Eq |A2l and Eq |A3l 
to 
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= (A18) 



and the unique solution 

Xq = Xi = -X 2 = -X3. 

Alternatively we can have: (2.ii) 
4tt 



(A19) 



(A20) 



z + z 3 = -- 



Zl + z 2 



> 2/o + y3 = = yi+y2 
(A21) 

which reduces the original equations Eq lAll and Eq |A3l 
to 
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= (A22) 



cos ( 2/2 - y J + cos \y° + y 

and the unique solution 

yo = -yi =yi = -2/3 



= (A23) 



(A24) 
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